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ABSTRACT 

We calculate the velocity and thickness of a bubble wall at the electroweak phase 
transition in the Minimal Standard Model. We model the wall with semiclassical 
equations of motion and show that friction arises from the deviation of massive particle 
populations from thermal equilibrium. We treat these with Boltzmann equations 
in a fluid approximation in the background of the wall. Our analysis improves on 
the previous work by using the two loop effective potential, accounting for particle 
transport, and determining the wall thickness dynamically. We find that the wall is 
significantly thicker than at phase equilibrium, and that the velocity is fairly high, 
v w ~ 0.7c, and quite weakly dependent on the Higgs mass. 



1 INTRODUCTION 

There has been a growing interest in the idea that the baryon asymmetry of the 
Universe may be created at a first order electroweak phase transition. However, 
ingredients needed to construct the whole picture are still missing. The exact nature 
of the Higgs mechanism is unknown, and the simplest model, the Minimal Standard 
Model, apparently does not contain sufficient CP violation for baryogenesis; we must 
consider extensions, such as the two Higgs model. Even in the Minimal Standard 
Model we have difficulties computing the finite temperature effective potential (which 
is needed to determine the strength of the phase transition) and the dynamics of the 
transition. 

Recently there have been advances in calculating the effective potential; the two 
loop contribution has been evaluated |], ||, and there has been progress in under- 
standing nonperturbative effects ]3|, ||, Both results support the view that the 
phase transition is first order and strong enough to proceed by bubble nucleation and 
growth (even when the Higgs mass is moderately large, win ~ raw)- 

To model baryogenesis accurately one also needs to know the profile and velocity 
of an expanding bubble wall. The wall velocity is friction limited; but determining 
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the strength of frictive effects involves determining the non-equilibrium populations 
of massive particles in the vicinity of the wall, which is difficult || [7|, |8| . The recent 
discovery that the top quark is very heavy |J suggests that top quarks may be the 
dominant source of friction, in which case an approximation which models top quarks 
with good accuracy should improve our understanding of the wall motion. In this 
letter we re-analyze the bubble wall's velocity and its shape, working in the Minimal 
Standard Model and using the fluid approximation to model the particle populations, 
which should treat fermions fairly well. The technique also allows us to account for 
transport and to determine the wall thickness dynamically (taking a specific Ansatz 
for the wall profile). 



2 Equation of Motion for the Higgs vev 

We intend to study the dynamics of infrared condensates in the Higgs field $. Such 
condensates should behave semiclassically to a good approximation. From the terms 
in the Electroweak Lagrangian containing <3>, 



(where the sum runs over the quarks and leptons if), and y denotes the Yukawa 
coupling), we find the equations of motion for <p (where & = [0, 0/ a/2]) to be 

2 

□0 + V\4>) - ^TtA 2 + ig w A»d,<P + + £ -tufa = (1) 

Here A and if> are quantum operators. (For simplicity we have set the Weinberg angle 
tan#ty = 0.) It is reasonable to take thermal averages of these operators using WKB 
wave functions. This is because the wall will be much thicker than the thermal length 
T" 1 , which characterizes the reciprocal momenta of particles in the plasma. In this 
approximation we get 

□, + m + E £/_£L_ /(ill) = (2) 

where V is the renormalized vacuum potential, / is the phase space population density 
(in the background of a propagating wall) and the sum includes all massive physical 
degrees of freedom. Note the condensed notation: m = y<f>J \/2 for quarks and leptons 
and m = g w <p/2 for the gauge fields. 

To model the population density / we assume a small departure from the equilib- 
rium population / and write f = fo + Sf. The vacuum contribution V'(<p) and the 
contribution from fo combine to give the finite temperature effective potential Vp(<f>). 
Thus we have 



* + + E^/p$b"(»*> = o (3) 

We see that the frictive force arises due to the departure from thermal equilibrium 
5f. We will use this equation, an expression for Vp, and equations for Sf to compute 
the wall velocity and its shape once it has reached a planar steady state, = <p(z + 
v w t),8f = 8f(z + v w t). 
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3 Effective Potential 



The high temperature expansion of the one loop effective potential is |7| 



V T {(j>) = D(T 2 - T 2 )0 2 - E(f> 3 T + -^0 4 



(4) 



with D = (2m 2 + 2m 2 w + m|)/8w 2 ~ .167, E = (2m 3 w + m 3 z ) / (4ttv 3 ) ~ 0.01, A T the 
Higgs self coupling at a scale roughly given by T, and T = [m 2 H + 3(4m 4 — 2m 4 y — 
m|)/87r 2 fQ]/4D. Recently the authors of [|TJ and 0] have computed the two loop 
expression. The most important changes are that E becomes (4m^ + 2m 3 z + 3(1 + 
"\/3)A 3 / 2 ) / (127T1>q) and a qualitatively new term —Bcf) 2 T 2 log(0/T) appears. The result 
of is5~ (1.4^ -.48^v / A + 5.l5f 2 A-6.7A 2 )/(167r 2 ). It may also be important to 
include nonperturbative effects. Shaposhnikov proposes a term — (A fg^T 4 / 12) Pit(0) 
0. The function Pit describes the contribution of a gauge condensate and it is roughly 
constant near = and falls exponentially as exp —<f)/g 2 T for large (f>. We add such 
a term to parametrize our ignorance about the free energy of the symmetric phase; 
we use — (A fg^T 4 / 12) sech(320/T) because it is simple - the exact form of Pit will 
have little effect on our calculations - and we take At and Af as unknowns. 

What should we use for T? Most of space is converted to the broken phase by 
bubbles which nucleate when the critical bubble free energy reaches S ~ 100T [|TTJ. 
We have computed S using our form for V T and the techniques used in [|, 0, || and 
find, for instance, that for A T = 0.03, A f = 0, and S = 100T that T = 1.00006T . 
(Note that because of the B term, T is no longer the spinodal temperature.) 

4 Fluid Equations 

Next we must determine 5f, the deviation from equilibrium in the presence of the 
moving wall. Our starting point is the Boltzmann equation, 



v z = Pz/E is the velocity (with the broken phase at positive z), and p z = —d z E is the 
force on the particle. In a complete description each particle species in the plasma 
would be described with a Boltzmann equation. We allow ourselves the approximation 
that all species but top quarks and perhaps W bosons are in equilibrium and we 
neglect the slight change in the background temperature across the wall. The former 
approximation is reasonable as the induced deviation from equilibrium goes as m 2 and 
top quarks and W bosons are the heaviest particles. We can correct for the latter 
approximation by using the work of Enqvist et al. |12| and Heckler jnj, who have 
used hydrodynamic conservation laws to relate wall velocity to entropy production; 
for simplicity we will not do so here. We also use a single / for tops and antitops of 
both helicity. In the Minimal Standard Model this is reasonable as there is almost no 
CP violation, and the difference in transport properties arises only at the subleading 
level of weak scatterings. In two doublet models we might need to be more careful. 

The Boltzmann equations are nonlinear partial integro-differential relations and as 
such are analytically intractable. To solve them we make the fluid Ansatz, assuming 
/ to be of the form 




I 2 is the particle energy, 



(5) 
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f , E-E8T/T-p z v-u (6) 
1 + CXp hjr-* 1 

where we have written explicitly three types of perturbations: chemical potential ^, 
temperature ST and velocity v. This Ansatz is a truncation of an expansion in powers 
of momentum; it gives a reasonable description of the populations of thermal energy 
particles when the background varies slowly on the scale of the diffusion length. For 
top quarks this should be sufficient as the diffusion length is short and the influence 
of infrared particles is phase space suppressed. For W bosons, Bose statistics give 
large infrared particle populations, and the fluid approximation is unreliable unless 
W bosons thermalize on time scales short compared to their annihilation rate. We 
consider the fluid equations for W bosons as a guide for their importance and con- 
centrate on top quarks. We will furthermore work to linear order in perturbations 
which are of order (m/irT) 2 , and therefore naturally small. 

With a three parameter Ansatz (|6p we cannot ask that the full Boltzmann equa- 
tions be satisfied, but only impose that three moments be satisfied, namely the in- 

3 3 3 

tegrals over f (t p, JclpE, and J d p p z . Working in the fluid frame and using 



d t f(z + v w t) = v w f we obtain the following equations |T4 



^ ST' 1 . ^ _ a _ ST 

av w ^+v w — + -v' + F 1 = -T^-T T1 — (7) 

bv^+vJ-^ + h' + F, = -T^-vJ-^ (8) 

b^ + — + v w v' + = -T v v (9) 

where a = 2( 2 /%3, b = 3n /4p = 3^3/14(^4, C is the Riemann (^-function: ( 2 = k 2 /6, 
Cs w 1-202, C 4 = vr 4 /90, n = 3C 3 T 3 /4vr 2 , p = 21^/^. 
The force terms are: 

v w In 2 (m 2 )' v w ( 2 (m 2 )' 

1 9C 3 T 2 ' 2 42( 4 T 2 1 } 

They drive the plasma out of equilibrium, while scatterings restore it. Scatterings 
keep the plasma near equilibrium when the wall passage time L/v w ^> T -1 . 

We have computed the coefficients on the r.h.s. of (0) - (|D including all diagrams 
which contribute to order ot 2 \oga s . We found |L5| Y^i m T/26, Tti ~ T/14, T^ 2 ps 



T/60, T T2 ps T/16, r„i ~ T/13, where we used a s = a s (m z ) — 0.12. For W bosons 
the values are about half as large. 

The force terms are proportional to m 2 , so in general Sf oc m 2 . Note that the 
friction term in Eq. (^) is proportional to m 2 Sf cx m 4 . For bosons the coefficients 
in Eqs. (0) - (|10D differ, and in particular F% oc m 2 log(2/m). The integral in the 
equation of motion involving fi also contains a log enhancement; but because mf/m^ 
is very large, W bosons still produce less friction than top quarks. 
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5 Computing Velocity and Profile 



With our Ansatz for / we can rewrite Eq. (^) for top quarks (with —> 12) as 

>ln2 + C 2 £T" 



=o (n) 



4tt 2 

This equation and the fluid equations form a system of nonlinear differential equations 



for the wall profile and velocity. We will attempt to solve them in ||15|| , but here we 
will content ourselves with an Ansatz for 6, 



<P=f (^1 + tanh^-^j (12) 

where 0o is the value of <fi in the asymmetric phase and v w the wall velocity and L 
the wall thickness in the plasma frame are treated as undetermined parameters. This 
Ansatz is chosen because the static equilibrium wall shape in the one loop approxi- 
mation is of this form. 

Again, because we have restricted the form of we cannot ask that the full 
equations of motion be satisfied; we can only enforce two moments. The natural 
choices are the space integral of the equation of motion times d<p/dv w and d<p/dL. 
Note that d<p/dv w = t<p' and d<p/dL = —(x/L)(f)', so an equivalent set of conditions is 

J [Eq. {UMdz = , J [Eq. ^[jtffdz = (13) 

These equations have a simple physical interpretation. The first equation is the total 
pressure on the wall in its rest frame ||; if it were non-zero the wall would accelerate, 
changing v w . The second equation is the asymmetry in the pressure between the 
front and back edges of the wall; if it were nonzero the wall would be compressed or 
stretched, changing L. 

The integrals for the first two terms in ([TT]) are 

;n0 + ^(0))0' = v T (<f>o)-V T (p) = -AV T (14) 
□0 + ^(0)]-^/ = il_|M-I[AV r + H] (15) 
E = - 1)2- + - + (2,2 + \ in ^ (16) 

Note that the 0<p term acts to stretch the wall (increase L) while Vt acts to accelerate 
and compress the wall. The coefficient 2.82 in the last term is the only place where 
our choice for the function Pit enters our computation. 

We will first get a rough estimate of the wall velocity and thickness by solving Eqs. 
(0) ~ © an d dUD ignoring transport, by which we mean we will ignore the derivative 
terms on the l.h.s. of the fluid equations. Transport reduces friction because particles 
tend to flow off the wall, where they contribute less to equations (|T3|). We will also 
ignore 5T, which turns out to be a good approximation. The expression for /x now 
becomes rather simple: 

to 2 2 #' 

« = v -^ v 'T^f (17) 



Z + V w t 
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Note that /i does not depend on other decay rates apart from r M i. 

Using )(00') 2 = 0o/ loL , Jx(H') 2 = 0o/ 24 and © " & one can solve for 
L w = L/(l - vl) 1/2 and v w : 



1 1 fAV T n \ 157r 2 C 3 „ _ AVj 



T 



For m t = 174GeV, A T = 0.03 and A f = these give L w ~ 29T" 1 and 7^ ~ 1.1, a 
mildly relativistic and fairly thick wall. Note that L ~ 20T -1 is much thicker than 
the top quark diffusion constant D ~ 4T" 1 , so the fluid approximation is in good 
shape. 

Now we solve the problem including transport. First we must find the contribu- 
tions to involving \x and ST. This is most easily accomplished by Fourier analysis. 
Let us write Eqs. (0) - @ in a matrix notation, 

AS' + T5 = F<M>' (19) 

where 8 is a column vector of /x, <5T, t>, A is the matrix of coefficients for the derivative 
terms, T is a matrix of the decay constants, and F is a column vector of the coefficients 
for the force terms. Note that A is velocity dependent, and F is linear in velocity. In 
Fourier space fll9D becomes 

ik8 + A- 1 T8 = A- 1 F^)' (20) 

which may be solved by eigenvalue methods. Denoting the eigenvalues and eigenvec- 
tors of A _1 T as Xi and & and expanding A _1 F = a^, we find 

& = «8/2) (l - f ) (*§=) csch (22) 

This gives an explicit expression for p, and ST. 

Finally, we can convert the relevant integrals of the friction term 
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\pi(x) In 2 + C2ST '(x)]<f></)' \x)dx 



3yfT r r x 



r2 



,, , \ji(x)\n2 + C 2 ST(x)}-4)(f) l (x)dx (23) 

7T J L 

into k-space_ integrals using the relations / fi(x)f>i(x)dx = J fi(—k)f2(k)dk/2ii and 
xf(x) =r> idf(k)/dk. This yields integrals of form 

1 / k 2 L 2 \ fkLn , kLn\ 2 dk 



csch — (24) 

X + ik \ A J \ 2 2 J 2vr y J 

which may be converted to a rapidly converging infinite sum by residue integration, 
or performed numerically. 

This completes the evaluation of all terms in Eqs. (|13|). These equations each 
define a curve in the space of v w and L. The intersection of these curves is a self con- 
sistent solution for the wall shape and velocity within the Ansdtze and approximations 
we have made. 
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6 Results 



We have solved these simultaneous conditions for some representative values of Xt and 
Af. We find that the friction from W bosons, calculated in the fluid approximation, 
is about half that from top quarks. Though our techniques are different than those of 
||, we get a similar numerical value for the friction from Ws. We have included them 
in our analysis. Using the temperature where the critical bubble action S = 100T, 
we find 
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The last two columns are the velocity and thickness of the wall when we treat v w 
as a free parameter but fix L to the value derived from fll5"D without the contribution 
from friction. This value of L approximately equals the thickness at phase equilibrium. 
We see that the wall is significantly deformed by the frictive effects, and that this 
increases its velocity. Of course, if the deformation is large then we have little reason 
to believe that our wall shape Ansatz is accurate - one should model the shape more 
carefully than we have done here. The conclusion that the wall is fast and thick 
should be reliable, however. 

When we include a sizeable value of Af, the parameter describing nonperturbative 
symmetric phase effects, we find no solution. The two equations (|T3| ) turn out to be 
incompatible; the wall runs away, but maintains finite plasma frame thickness. This 
result probably comes from neglecting friction from the gauge condensate responsible 
for Af, which would compress the wall and prevent runaway. To remedy this short- 
coming we need a model for the nonequilibrium dynamics of nonperturbative infrared 
condensates. 

The situation in two doublet models may be quite different from what we have 
found here. In these theories there are several new massive (Higgs) bosons. The ones 
which do not couple to the top quark have quite long half-lives and sizable diffusion 
constants, and may be a major source of friction. 
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